Abstract. We report observations of the low-luminosity 
z = 5.50 quasar RD J030117+002025 (RD0301 hereafter) 
at 250 GHz (1.20 mm) using the Max-Planck Millimeter 
Bolometer (MAMBO) array at the IRAM 30-meter tele- 
scope. The quasar was detected with a 1.2 mm flux den- 
sity of 0.87 ± 0.20 mJy. The lack of detectable 1.4 GHz 
radio emission indicates that the millimeter emission is of 
thermal nature, making RD0301 the most distant dust- 
emission source known. When matching a 50 K grey body 
thermal far-infrared (FIR) spectrum to the observed mil- 
. limeter flux we imply a FIR luminosity « 4 x 10^^ L©, 
' which is comparable to the quasar's optical luminosity, 
i: If the FIR luminosity arises from massive star formation, 
■ the implied star formation rate would be ~ 600MQyr~^, 



comparable to that of the starburst galaxies which domi- 
nate the average star formation and FIR emission in the 
Qv^ ,; early Universe. The FIR luminosity of RD0301 is close to 
^—i ' the average of that found in optically far more luminous 
high-redshift quasars. The comparably high millimeter to 
optical brightness ratio of RD0301 is further evidence for 
that there is no strong correlation between the optical 
and millimeter brightness of high-redshift quasars, sup- 
■ porting the idea that in high-redshift quasars the dust is 
T— I ■ not heated by the AGN, but by starbursts. 
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1. Introduction 

Many higli-rcdshift quasars (QSOs) have recently been 
found through wide field imaging surveys, notably the 



ties or velocity dispersions in local galaxies indicates that 



Sloan Digital Sky Survey (SDSS; Schneider et al. 2001 



Anderson et al. 2001) and the Digital Palomar Sky Sur- 
vey (DPSS; Kennefick et al. |l9954|l995lj ; Djorgovski et 
al. |1999D . Now over 250 QSOs are known with redshifts 



2000a 



z > 3.6, thirteen of which are at z > 5 (Fan et al. 1999, 
2000b| ; Zheng e t al. pOOOj ; Stern et al. |2000| ; Sharp 



et al . 2001; Fan et al. 2001) and one at z > 6 (Fan et al. 
2001). All but one of the z > 5 QSOs were found from 
magnitude i « 21.5 limited surveys, limiting their implied 
rest frame blue magnitudes to M-q < — 26. 

In a 74 arcmin^ field observed for the SPICES sur- 
vey (Stern et al. |200lD , Stern et 
the redshift 5.50 QSO RD (i.e. 



al. ( pOOOD discovered 
"R-band Dropout" ) 
J030117-h002025 (RD0301 hereafter), which has an i-band 
AB magnitude of 23.4 and Mb ~ —23.4 in an Einstein-de 
Sitter Universe {Hq = 50 km s~^Mpc~^, go = 0.5), and 
Mb « -24.0 in a A-cosmology [Hq = 65 km s"^ Mpc~\ 
= 0.7, f2ni = 0.3) which we adopt throughout this pa- 
per. RD0301 is significantly fainter than the typical QSO 
found in the wide field SDSS and DPSS surveys. The dis- 
covery of such a very high redshift QSO was surprising 
because the surface density of such objects, as implied 
by the intermediate-redshift QSO luminosity function and 
redshift evolution models, is so low that the probability 
to find one in the observed field is only ~ 15% (Stern et 
al. 2000 ). The discovery may therefore imply that the lu- 
minosity evolution of faint high-redshift QSOs is weaker 
than expected, and that there could be a larger, optically 
faint QSO population at high redshifts. This possibility is 
particularly interesting in terms of which sources re-ionize 
the early Universe (cf., the recent Gunn-Peterson trough 
results of Becker et al. [ 2001 | and Djorgovski et al. | 2001 |). 

Another important question in studies of the early Uni- 
verse is how the growth of massive black holes relates to 
the formation of early stellar populations. A tight cor- 
relation between black hole masses and bulge luminosi- 



their formation was closely related (Magorrian et al. 1998). 
The recent detection of strong thermal dust emission from 
many high-redshift QSOs shows that vigorous star for- 
mation is coeval with black hole growth (Carilli et al. 
p001a| ; Omont et al. |2001| ; Isaak et al. pOOl] ). Interestingly, 
the FIR luminosity of these QSOs does not seem to be 
well correlated with their optical luminosity in the stud- 
ied range of optical luminosities, A/b = —26 to —29.5. A 
similar result was previously found also in lower redshift. 



lower optic al luminosity QSOs (Sanders et al. 1989; Chini 



et al. 1989). Such a lack of correlation might suggest that 



the warm dust is not heated by the QSO. Instead, star 
formation and black hole accretion could have had a com- 
mon cause, which might be the infall of material toward 
the center of the host galaxy. 

The QSOs observed in the millimeter and submillime- 
ter surveys are the optically most luminous objects, and 
their implied optical luminosities are on average ten times 
higher than their FIR luminosities, so that dust heating 
by the AGN cannot be excluded on energetic grounds. It 
is of great interest therefore to investigate whether strong 
millimeter emission is also found in high-redshift QSOs 
which are optically much fainter than those observed so 
far. In this respect, RD0301 is a unique object: it is the 
optically least luminous QSO known at z > 4, and it is 
one of the highest redshift quasars ever detected. 

Searching for high-redshift dust emission is of great 
interest also because for them there was little time to pro- 
duce the observed dust. At redshift 5.5, the universe was 
about 1 Gyr old, a time not much longer than the dynami- 
cal timescale of typical galaxies. Observing large amounts 
of dust at such redshifts thus sets back the first epochs 
of vigorous star formation irrespective of whether the ob- 
served dust was heated by an AGN or a starburst. 

2. Observations and results 

Our observations of RD0301 were made during the winters 
of 2000 and 2001 with the 37-channel Max-Planck Mil- 



Send offprint requests to: F. Bertoldi, bertoIdi@mpifr.bonn- 
mpg.de 



limete r Bolometer array (MAMBO; Kreysa et al. 1999a , 
1999b) at the IRAM 30- meter telescope on Pico Veleta, 
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Fig. 1. Time-averaged, opacity corrected, and sky-noise 
subtracted flux densities for 33 of the 37 MAMBO bofome- 
ters. Four disfunctional bofometers are ignored. The thick 
hne shows the time-averaged signal from the central, on- 
source bolometer, the grey-shaded envelope is reflects its 
Icr error. The thin solid lines are the average on-sky 
bolometer signals, and the dashed lines show the la scat- 
ter of the average on-sky bolometer signals. 



Spain. MAMBO is sensitive between 190 and 315 GHz, 
with half-power sensitivity limits at 210 and 290 GHz, 
and an effective bandwidth center for steep thermal spec- 
tra at ~ 250 GHz. The bolometer feed horns are matched 
to the telescope FHWM beam at 1.2 mm of 10'.'7. They are 
arranged in an hexagonal pattern with a beam separation 
of 22". The source was observed with the array's central 
channel, using the "on-off" observing mode, in which the 
telescope secondary mirror chops in azimuth by 30 to 50" 
at a rate of 2 Hz. For on-off observations the target is po- 
sitioned on the array's central bolometer, and after every 
10 seconds of integration the telescope is nodded so that 
the previous "off" beam, which typically misses the other 
bolometer beams, becomes the "on" beam. Thus for only 
half the on-sky observing time the target is positioned on 
any bolometer. A scan usually lasts for 12 or 16 such 10 
sec sub-scans, i.e., 6 to 8 on-off cycles. 

The telescope pointing was monitored regularly on 
two ~ 10 deg distant pointing sources, 0336—019 and 
02214-067, and was in most cases stable within 2". The 
signals from the pointing sources were monitored for sensi- 
tivity problems. The gain calibration was performed using 
observations of Mars, Uranus, and Ceres. We here adopt 
a calibration factor of 12 500 counts per Jansky for 0.5 
sec integration bins, a value which is expected to be ac- 
curate to within 20%. The point-source sensitivity was on 
average '~ 30 mJy s^/^, i.e., in one second of on-sky (on 



and off source) integration, an rms noise level of 30 mJy 
is reached. 

The data were analyzed using the MO PSI so ftware 
package created by Robert Zylka (see Zylka 1998| ). Cor- 
related sky noise was subtracted from each channnel as 
the weighted mean of the signals from the best-correlating 
surrounding channels. 

RD0301 was observed on seventeen different nights be- 
tween February 2000 and March 2001. The atmospheric 
conditions were mostly very good, with zenith opacities 
between 0.1 and 0.2. From a total of 184 scans, each with 
120 or 180 sec on-sky exposure, we eliminated one obser- 
vation (1 March 2000) of 15 scans because of a 5" pointing 
drift. Six other scans were eliminated throughout because 
they showed unusual noise, possibly due to passing clouds. 
To improve the quality of the data we furthermore ig- 
nored 10 scans at elevations below 27 degrees, where the 
gain-elevation and opacity corrections become large and 
uncertain. 

The total on-sky integration time was 6.0 hours. The 
scans were reduced individually, and the signals were then 
averaged, which yields a target flux of 0.87 ± 0.20 mJy 
(Fig. I). If the scans with elevation < 27 deg were in- 
cluded, the signal rises to 0.93 ± 0.21 mJy. When weight- 
ing the scan signals by each scans' inverse average mean 
square noise level, the signal becomes 0.76±0.19 mJy. The 
dispersion in the integrated signals of the off-target chan- 
nels is 0.26 mJy. The source signal appears stronger in 
the winter 1999/2000 observations, which might in part 
be due to the fact that many of these observations were 
performed at low elevation, making them subject to larger 
gain-elevation and opacity uncertainties. 

The coadded data of both years show an average signal 
that stabilizes at a significant level compared with the tar- 
get channel's integrated noise, and also compared with the 
signal of the off-target channels. We therefore believe that 
the source is detected with the confidence level suggested 
by the errors quoted. 

To our knowledge this is the deepest integration ever at 
millimeter or submillimeter wavelengths, and the faintest 
thermal emission ever detected. As a promising outlook 
we note that this object would be detected by ALMA in 
a few seconds, with an angular resolution of less than 0.1 
arcsec. 



3. Discussion 

3.1. Thermal nature of the millimeter emission 

At redshift 5.5 our observing frequancy of 250 GHz cor- 
responds to an emitted frequency of 1625 GHz or a wave- 
length of 185 /im. From a flux measurement at only one 
frequency we cannot distinguish whether the emission fol- 
lows a steep thermal spectrum, or constitutes the millime- 
ter extension of a radio synchrotron spectrum. The VLA 
1.4 GHz FIRST survey (Becker, White & Helfand |l995|) 
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however yields a 5(7 upper limit flux density of c:± 1 mJy 
within 30" from RD0301. Since rising non-thermal spectra 
are not observed at high radio frequencies, it is unlikely 
that the millimeter emission is synchrotron radiation. Car- 



illi et al. (2001a, 20011:) show that for most of the high- 
redshift QSOs detected at 1.2 mm warm dust is the likely 
source of the emission. 



3. 2. FIR luminosity, dust mass, and star formation rate 

We therefore conclude that the millimeter flux of RD0301 
is thermal in nature. When matching the observed 1.2 mm 
flux of 0.87 mJy to a grey body with a dust temperature 
of 50 K and an emissivity index /3 = 1.5 (well adapted 
for high-redshift sources - see, e.g., Benford et al. 1999), 
the FIR luminosity of RD0301 is ipiR = 4.0 x IO^^l©, 
and the implied dust mass k, IO'^Mq. The FIR lumi- 
nosity of RD0301 is comparable to its optical luminosity, 
Lo-pt ~ 3.7 X 10^^ L0, given by the rest frame blue magni- 
tude Mb w —24.0 (for a A-cosmology). We here adopted 
the blue luminosity bolometric correction factor of ^ 12 
derived for the PG sample by Elvis et al. ( 1994 ). 

Both the FIR luminosity and the dust mass are compa- 
rable to those derived for local infrared-luminous galaxies 
such as Arp 220 and Mrk 231 (e.g.. 



Radford et al. |1991D . 
The thermal luminosity of RD0301 is comparable also to 
that of the bright high-redshift starburst galaxies found in 



deep optical surveys (Adelberger & Steidel 2000), and of 
those thought to be producing the bulk of the FIR back- 
ground emission (Peacock et al. pOOOj ). If the FIR luminos- 
ity of RD0301 arose from a continuous starburst of age 10 
to 10 Myr with a modified Salpeter IMF (see Omont 
et al. 2001 for details on the adopted starburst model), 
the implied star formation rate would be ~ 600 M© yr~^. 
This value again is similar to those inferred for Arp 200 
and Mrk 231. 



3.3. Comparison with the SDSS and DPSS QSOs 

The millimeter flux of RD0301 is typical for the optically 
brightest high-redshift QSOs. When averaging (as if we 
observed one single object) the MAMBO observations of 
the 112 re dshift 3.6 to 5.0 QSOs obs erved in the Carilli et 
al. (^OOlaD SDSS and Omont et al. (|200l|) DPSS surveys, 
one obtains an average flux density of 1.3±0.05 mJy, which 
implies an average luminosity of LpiR ~ 5 x 10^^ Lq. The 
dust luminosity of RD0301 is therefore close to the average 
dust luminosity of a sample of quasars with — 26 > ilfe > 
—29.5, sources which are optically 10 to 100 times brighter 

than RD0301. 

The mean blue magnitudes of the Carilli et al. ( ^OOla ) 
and Omont et al. (1 200 1| ) QSOs are -27.0 and -27.9 mag, 
respectively (Fig. |). Although the QSOs of these sam- 
ples span a range of about 25 in optic al bri ghtness, the 
DPSS QSOs studied by Omont et al. (^00 1[ ) are on av- 
erage twice as bright in the optical as the SDSS QSOs 
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Fig. 2. MAMBO 1.2 mm flux density as a function of blue 
magnitude for the z > 3.6 QSOs observed in the Carilli 
et al. (open symbols) and Omont et al. (filled symbols) 
surveys. Average source fluxes for the respective samples 
are indicated as large open squares, which are placed at 
the median Mb of the respective sample. 



observed by Carilli et al. ( 2001a ). Similarly, the average 
1.2 mm flux of the QSOs in the respective samples differs 
by a factor two, 0.92 mJy versus 2.1 mJy (large squares 
in Fig. H). However, the significance of a possible corre- 
lation between the average optical and millimeter bright- 
ness of high z QSOs remains low since the scatter in both 
quantities is much larger than the difference between the 
samples' averages (see Omont et al. |2002] for a detailed 
statistical discussion) . 

RD0301 does not follow the possible trend of increas- 
ing millimeter flux with optical luminosity seen in the 
QSO millimeter surveys, as Fig. || well illustrates. Fur- 
ther (sub) millimeter observations of optically faint {Mb > 
—26) high-redshift QSOs are needed to clarify the relation 
between their thermal and optical luminosities. A lack of 
such a correlation would support the idea that the dust 
heating is caused by young, massive stars, and not by the 
AGN. This hypothesis is also supported by the detection 
of CO emission in a number of high-redshift QSOs (Guil- 
Cox et al. 



loteau 2001 



2002), and by the lack of strong 



radio emission from most of them (Carilli et al. ^OOlb ) 



Alternatively, unification models of AGN might not 
necessarily predict a strong correlation between the 
orientation-dependent optical emission and the isotropic 
FIR flux from the AGN. But still the AGNs' FIR flux 
should correlate well with its hard X-ray emission, which 
does not suffer much from extinction. Millimeter-studies 
of high-redshift, hard X-ray sources would be most usefu l 
to address this interesting issue (e.g. Page et al. 2001). 
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SIRTF observations of mid-IR PAH features should also 
be able to clarify the origin of the FIR emission: a strong 
correlation between PAH and FIR emission supports the 
starburst origin of the FIR luminosity (Genzel et al. 19981 ). 
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